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Abstract 



Wc study uniaxial energy minimizcrs within the Landau-de Gennes theory for nematic liquid 
crystals, subject to dirichlet boundary conditions. Topological defects in such minimizcrs cor- 
respond to the zeros of the corresponding equilibrium field. Wc consider two-dimensional and 
three-dimensional domains separately and study the correspondence between Landau-de Gennes 
theory and Ginzburg-Landau theory for superconductors. We obtain results for the location and 
dimensionality of the defect set, the minimizer profile near the defect set and study the qualitative 
properties of uniaxial energy minimizers away from the defect set, in the physically relevant case 
of vanishing elastic constant. In the three-dimensional case, we establish the C 1, "-convergence of 
uniaxial minimizers to a limiting harmonic map, away from the defect set, for some < a < 1. 
Some generalizations for biaxial minimizers are also discussed. This work is motivated by the 
study of defects in liquid crystalline systems and their applications. 

1 Introduction 

Nematic liquid crystals are examples of mesophases whose physical properties are intermediate 
between those of a typical liquid and a crystalline solid [10] . The constituent rod- like molecules have 
no translational order but exhibit a certain degree of long-range orientational ordering. Consequently, 
liquid crystals are anisotropic media and this makes them suitable for a wide range of physical 
applications and the subject of very interesting mathematical modelling [13J . 

The Landau-de Gennes theory is a general continuum theory for nematic liquid crystals |10[ [28] . 
It describes the state of a nematic liquid crystal by a symmetric, traceless 3x3 matrix - the Q-tensor 
order parameter, that is denned in terms of anisotropic macroscopic quantities, such as the magnetic 
susceptibility and the dielectric anisotropy. Nematic liquid crystals are said to be in the (a) biaxial 
phase when Q has three distinct eigenvalues, (b) uniaxial phase when Q has a pair of equal non-zero 
eigenvalues and (c) isotropic phase when Q has three equal eigenvalues or equivalently when Q = 0. 
For a general biaxial phase, Q can be written in the form |15[ [20] 
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where s, r are scalar order parameters, n, m are eigenvectors of Q and I is the 3x3 identity matrix. 
In the uniaxial phase, Q takes the simpler form of 



where s is a scalar order parameter that measures the degree of orientational ordering about the 
distinguished eigenvector n. 

The Landau-de Gennes energy functional, leg, is a nonlinear integral functional of Q and its 
spatial derivatives. In the absence of any surface energies or external fields, leg is given by JTOJ |2"U] 



where f2 is the domain, /b(Q) is the bulk energy density that dictates the preferred phase - isotropic, 
uniaxial or biaxial, L is a positive material-dependent elastic constant and |VQ| 2 is an elastic energy 
density that penalizes spatial inhomogeneities. The equilibrium, physically observable configurations 
correspond to either global or local Landau-de Gennes energy minimizers, subject to the imposed 
boundary conditions. 

In this paper, we study uniaxial global minimizers (of the form @) of the Landau-de Gennes 
energy functional. There is the very important underlying question - do uniaxial global minimizers 
actually exist [15] ? This is an open question but there is substantial numerical and experimental 
evidence to show that global energy minimizers are largely uniaxial almost everywhere, in the sense 
that they have a small degree of biaxiality. Therefore, a rigorous study of uniaxial global minimizers is 
the first step in the mathematical analysis of arbitrary global minimizers and the interplay between 
biaxiality and uniaxiality. Secondly, stable uniaxial configurations do exist, at least for certain 
temperature regimes and certain physically realistic choices of the material-dependent constants 
P3S E31 [25] . Some of our results extend to stable uniaxial configurations with appropriately bounded 
energy and can then be used to understand these configuration structures and the nature of their 
singularities. We also point out that there are two widely-used continuum theories for purely uniaxial 
liquid crystal phases - the Oseen-Frank theory and the Ericksen theory, both of which have received 
considerable attention amongst the mathematical analysts and the numerical modellers O [T2| HH1 [8] . 
In fact, uniaxiality is one of the most frequently used assumptions in the theoretical study of liquid 
crystalline systems, even in the context of applications. 

We assume that uniaxial global minimizers exist for each L > and then establish various 
properties of such uniaxial global minimizers throughout the paper. The paper is organized as 
follows. In Section [2j we introduce some basic notation and terminology. In Section O we study 
Landau-de Gennes minimizers on two-dimensional (2D) domains and establish a 1 — 1 correspondence 
between Landau-de Gennes theory and Ginzburg-Landau theory. In Section [H we recall useful 
results from |15| that are crucial for the development of this paper. We study uniaxial Landau-de 
Gennes minimizers on three-dimensional (3D) domains in the low-temperature regime. There are 
important differences between the 2D and 3D cases and the standard Ginzburg-Landau techniques 
do not extend to the 31? case. We derive the governing equations for uniaxial global minimizers and 
obtain qualitative information about topological defects. The scalar order parameter V necessarily 
vanishes at the defect locations. The defect locations are prescribed in terms of the singular set of 
a limiting harmonic map and using asymptotic methods, we show that the leading eigenvector n 
(see ([2])) necessarily has a radial-hedgehog type of profile in the immediate neighbourhood of each 
isolated point defect. Our result is analogous to a powerful result on singularity profiles in [5], 
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where the authors work within the Oseen-Frank theory for uniaxial liquid crystals with constant 
order parameter s. In Section \5\ we study the qualitative properties of uniaxial global minimizers 
away from the defect set, in the limit L — > + . The elastic constant L to typically several orders 
of magnitude smaller than the other material-dependent constants and hence the L — > limit is 
physically realistic [22] . We adapt the small energy regularity theorem of [7| to the Landau-de Gennes 
framework and prove the C 1,Q -convergence of uniaxial global minimizers to a limiting harmonic map, 
away from the defect set, as L — > 0. This convergence result encodes quantitative information about 
the corresponding scalar order parameter. In Section [6j we discuss various generalizations of our 
results to uniaxial solutions with bounded energy and to the completely general biaxial case. The 
uniaxial case in 3D can be viewed as a generalized Ginzburg-Landau theory from K 3 to M 3 although 
there are important technical differences. However, the biaxial case presents a whole host of new 
mathematical difficulties; there are five degrees of freedom in the biaxial case and the additional 
degrees of freedom give us more possibilities, particularly in the context of defects. The methods 
in this paper contribute to the development of a generalized Ginzburg-Landau theory from M 3 to 
higher dimensions (M 5 in this case), for non-standard non-convex multi-well bulk potentials. 



2 Preliminaries 

Let 5 C M 3x3 denote the space of symmetric, traceless 3x3 matrices i.e. 

S d = f {Q G M 3x3 ; Qij = Q ju Q« = 0} 

where we have used the Einstein summation convention; the Einstein convention will be used in the 
rest of the paper. The corresponding matrix norm is defined to be 



|Q| v ^ r Q^=yQ-Q- z,j = 1...3. 

We take our domain £1 to be either a two-dimensional or three-dimensional bounded, connected 
and simply-connected set with smooth boundary, <9f2. We work with the simplest form of the bulk 
energy density, Jb, in ^ that allows for a first-order nematic- isotropic phase transition [2D]. We 
focus on the low-temperature regime; the function fs is bounded from below and can be written as 

MQ) = -ytr (Q 2 ) - b \ (Q 3 ) + | (tr(Q 2 )) 2 + C{a\ b 2 , c 2 ) (4) 

where a 2 ,b 2 ,c 2 G R + are material-dependent and temperature-dependent positive constants and 
C(a 2 , b 2 , c 2 ) is a positive constant that ensures /b(Q) > for all Q-tensors. We note that C(a 2 , b 2 , c 2 ) 
plays no role in energy minimization, in either spatially homogeneous or inhomogeneous cases. The 
negative coefficient of tr (Q 2 ) incorporates the fact that we are working below the nematic-isotropic 
transition temperature where Jb attains its minimum on the set of uniaxial Q-tensors given by [16] 



Qmm = |Q € S,Q = s + (^n®n- h^j X (5) 

with n G S 2 and 

b 2 + V& 4 + 24a 2 c 2 
S+ = 4^ ■ (6) 
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We study uniaxial global minimizers of leg , with fs as in (jH , with strong anchoring conditions 
or dirichlet boundary conditions. The prescribed boundary condition is given by 

Qb = s+ (n b ® n 6 - (7) 

where S W 1,2 {yi;M) (M = S 1 in 2D and M = S 2 in 3D) is a unit-vector field with non-zero 
topological degree d, when viewed as a map from dQ to M. Clearly, Q b G Q m in where Q m i n has 
been defined in ([5]). We define our admissible space to be 

A Q = {Qe W 1 ' 2 (Q; S);Q = Q b on dSl, with Q fe as in ©} , (8) 

where W 1 ' 2 (Q; S) is the Sobolev space of square-integrable Q-tensors with square-integrable first 
derivatives [§] . The existence of global energy minimizers for leg in the admissible space Aq follows 
readily from the direct methods in the calculus of variations |16|. 1 15j . For completeness, we recall that 

1 /2 

the W /1,2 -norm is given by ||Q||vk 1 . 2 (Q) = {Jq IQ| 2 + | VQ| 2 dx) . In addition to the VF 1,2 -norm, 
we also use the L°°-norm in this paper, defined to be ||Q||l°o(q) = ess sup xg0 JQ(x)| . 

Finally, we introduce the concept of a "limiting uniaxial harmonic map" Q° : VL — > Q m j n ; Q° is 
defined to be 

Q° = s+ (n 8) n - ±lj (9) 
where no is a minimizer of the Dirichlet energy 

I OF [n] = [ |Vn| 2 dV (10) 
Jn 

in the admissible space 

A n = {ne W 1 ' 2 (O; M) ; n = n b on dn} (11) 

where M = S 1 in 2D and M = S 2 in 3-D. The terminology limiting harmonic map stems from the 
fact that no is a harmonic unit-vector field [23] and it can be shown that Q° is a global minimizer of 
leg in the restricted class Aq n {Q m in} [SLPS]- We use the limiting harmonic map Q° to study the 
inter-relationship between the Landau-de Gennes theory and the Oseen-Frank theory for nematic 
liquid crystals. The Oseen-Frank theory is the simplest continuum theory for nematic liquid crystals, 
restricted to uniaxial phases with constant scalar order parameter [TO]. Working within the one- 
constant approximation, the Oseen-Frank energy reduces to the Dirichlet energy in ()10p and no, and 
hence Q°, is a global Oseen-Frank energy minimizer in the admissible space _4. n . 



3 The 2D case 

Let Q C K 2 be a bounded, connected and simply-connected domain with smooth boundary. Let S2 
denote the space of symmetric, traceless 2x2 matrices. Then Q £ S2 can be written as 

Q = A(n<g>n-m(g>m) (12) 

where A G R and n, m are the two orthonormal eigenvectors of Q. We note that there are only 
two degrees of freedom in the representation (|12p and hence we can think of Q as being a map 
Q : 17 — > M 2 . Using the identity, = n^n,- + mjnij, we can re-write (fl~2|) as 

Q = 2Afn®n-in (13) 
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where I is the 2x2 identity matrix. Thus, all admissible Q-tensors in two dimensions necessarily 
have a uniaxial structure as in §2§ 0. 

Straightforward calculations show that 

|Q| 2 = 2A 2 

trQ 3 = QijQjpQpi = i,j,p = 1,2. (14) 
Then the Landau-de Gennes energy functional in (|3j) simplifies to 

2«[Q] = /Jl v QI 2 + i{-y tr Q 2 + T( tr Q 2 ) 2 } dv 

for two-dimensional domains. The corresponding Euler-Lagrange equations are - 



(15) 



Q JJ -,fcfc = ^(-a 2 + c 2 |Q| 2 )Q ii i,j = 1,2 (16) 



and using the scaling Q = \/^a-Q, we obtain the following system of partial differential equations 

,2 



~ a 



Q 



2 



1 Qii i,j,k = 1,2 



Q = 2 ^n fe <g> n fe - ilj on 50. (17) 

This is identical to the Ginzburg-Landau equations for superconductors in two dimensions [26] and 
we are interested in the asymptotic properties of global energy minimizers either in the limit a 2 — > oo 
or L ->■ 0+. 

Let Q L be a global minimizer of leg m (|15p. in the admissible space 
„4 Q = {Q e W 1 ' 2 (0; 5 2 ) ; Q = s+ (n 6 <8> n b - |l) on 90} for a fixed L > 0. Then Q L is necessarily 
of the form 

Q L (x) = s L (x) ( n L (x) n L (x) - h) (18) 



for some scalar function s L : -> E and n L G W 1 ' 2 (fi; Let Q L = {x G 0; s L (x) = 0} denote 
the isotropic set of Q^. We have a topological^/ non-trivial boundary condition in (|17p . since nj, has 
non-zero topological degree when viewed as a map from dQ to S 1 . Hence, the unit- vector field n L 
necessarily has interior discontinuities and let S n denote the defect set of n^. Then 

s n c e L (19) 

and in what follows, we use existing results in the mathematical literature for Ginzburg-Landau 
theory in two dimensions, to make predictions about the structure and location of the isotropic set 
and the far-field properties of global energy minimizers. 

Dimension of @l [21 EQ: The isotropic set @l consists of |eZ| isolated points, {ai, . . -aui}, where 
d is the topological degree of the boundary condition in ([?])■ 



1 We can also think of Q £ 5*2 as being a symmetric, traceless 3x3 matrix: Q = (A+|)n®n+(| — A)m®m— |zcg>z G S 
where z is the unit-vector in the z-direction and S is the space of symmetric, traceless 3x3 matrices. 
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Defect locations [3 ], 126] : The configuration (a\, . . . aun minimizes the renormalized energy W 
over (bi, . . . bu\) £ Q^, which is defined by 

W(b u ...b M ) = -2TvJ2 lo s\bi-bj\-2TTj2 R ( b M ( 2 °) 

where i?(x,y) = \£(x, y) — log |x — y|,x, y £ M? and ^(x, y) is given by the solution of an explicit 
boundary-value problem. 

Far-field behaviour [3], [26]: Let {Q Lfc } denote a sequence of global energy minimizers for (|15p . 
where — > + as k — > oo. Then (up to a subsequence), 

Q Lfe -> Q* in C 1 ' (n \ G L ) , Va < 1 and in VF 1,p (ft), Vp £ [1, 2) 

for some Q* £ ni< p< 2VK 1,p (f^S 1 ). The limit Q* is the canonical harmonic map associated with 
a±, . . . a\ d \ and the degrees sgnd, . . . sgnd. 

The interested reader is referred to [3J [26] for the proofs. 

4 Uniaxial minimizers and their defect sets in 3D 

Let (] C I 3 be a bounded, connected and simply-connected domain with smooth boundary. An 
arbitrary Q-tensor field, Q : — > S can be written as 

3 

i=l i 

where e^ are the orthonormal eigenvectors, Aj are the corresponding eigenvalues and trQ 3 ^ in 
general. 

We study uniaxial global minimizers of the Landau-de Gennes energy functional, leg in ©, in 
the admissible space ^4q = {Q £ W 1,2 [Q; S) ; Q = Qf, on <9S7}. The corresponding Euler-Lagrange 
equations are 

LAQy = -a 2 Qi, - b 2 (Qi k Q kj - ^tr(Q 2 )^) + C 2 Q^tr(Q 2 ) i,j = 1,2,3, (21) 

where the term 6 2 ^-tr(Q 2 ) is a Lagrange multiplier associated with the tracelessness constraint. 
It follows from standard arguments in elliptic regularity that a global minimizer Q* is actually a 
classical solution of (|2ip and Q* is smooth and real analytic on 0, up to the boundary [15j . 

We assume that a uniaxial global minimizer exists for each L > 0. For Q uniaxial (of the form 
Q = s(n (g) n — |l) where s : O — > M. and n £ W 1,2 (Q; S 2 ), see ([2])), a direct calculation shows that 

(cUQkj - ^tr(Q 2 )) = ^Qij 
and hence, the Euler-Lagrange equations ()2ip simplify to 

LQij,kk = 3 (2c V - b 2 s - 3a 2 ) Q tJ , i, j = 1 . . . 3. (22) 
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Let Q L denote a uniaxial global Landau-de Gennes minimizer for a fixed L > 0. Then Q L is a 
classical solution of (|'22|) and we are interested in the qualitative properties of in the limit L — > 0. 

We briefly comment on limiting harmonic maps in a 3D setting: Q° = s+ (no <8> no — |l) where 
no is an energy minimizing harmonic map in the admissible space A n = {n € W 1,2 (Q; S 2 ); n = n& on <9f2} . 
Let Sq denote the singular set of no (and hence, of Q ). Then So consists of precisely \d\ isolated 
point singularities [5l [21] . 

We, next, quote important results from [15\ [T6] that are crucial for the analysis in this paper. 

Maximum principle [16] : Let Q be an arbitrary solution (not necessarily uniaxial) of the Euler- 
Lagrange equations (|2~T]) in the space Aq. Then 

HQlU°°(n) < \/f s + ( 23 ) 

where s+ has been defined in 

Strong convergence to Q° |15j : Let S! C I 3 be a bounded, connected and simply-connected 
domain with smooth boundary. Let {Q Lfe } be a sequence of uniaxial global minimizers of leg hi 
the admissible space .Aq ( leg and Aq have been defined in ([3]) and © respectively) where Lf, — > 
as k — > oo. Then Q Lk — > Q° strongly in W 1 ' 2 (O; S) (upto a subsequence), where Q° has been 
defined in ©. 

Interior and boundary monotonicity lemmas |15| : Let Q be an arbitrary solution of the Euler- 
Lagrange equations (|2ip . Define the normalized energy on balls -B(x, r) C O = {y G £1 : |x — y| < r}: 

HQ,*,r)= 1 -[ I|VQ| 2 + MQ) dV . 

r JB(^t) 2 L 

Then we have the following interior monotonicity lemma: 

J"(Q,x,r) < T{Q,x,R) Vx G U; r < R and B(x, R) C 0. (24) 

Similarly, for xo G 00, we define the region Sl r = Sin -B(xo,r) with r > 0, and the corresponding 
normalized energy to be 

£(Q,x ,r) = i / ^|VQ| 2 + M^dy. 
Then there exists ro > so that 

^-£> -C(a 2 ,b 2 ,c 2 ,Qb,r ^) < r < r (25) 

where the positive constant C is independent of L. 

The proofs of (|24p and (|25[) follow a standard pattern using the Pohozaev identity; complete de- 
tails can be found in [TS] . An immediate consequence of the strong convergence and the monotonicity 
lemmas is the following: 

Convergence of bulk energy density away from So [15] : Let {Q ife } be a sequence of global 
Landau-de Gennes energy minimizers in the admissible space .Aq, where — > as k — > oo. Assume 
that we have a sequence Q Lk with — > as k — > oo, such that Q ifc — > Q° in W 1 ' 2 ^, S), as k — > oo, 
where Q° has been defined in Q. 

For any compact set K C such that K contains no singularity of Q°, we have that 

lim / B (Q Lfc (x)) = xeK (26) 

L fe ^0 
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and the limit is uniform on K. 

Consider a sequence of uniaxial global Landau-de Gennes energy minimizers Q Lfc such that 
Lfc — > as k — > oo. Then Q ife can be written in the form 

Q Lk = s Lk (n Lk <g> n Lk - ^ij (27) 

for s k : Q — > M. and n k G W 1,2 (Q; S 2 ). Then (|26|) implies that (up to subsequence), s k converges 
uniformly to s+ everywhere away from Sq i.e. we have 

s fc (x)-s+| <e(L fc ,x) xefi\B i (S ) (28) 

where e — > + as k — > oo, 12,5 (So) is a small ^-neighbourhood of the singular set So and < 5 < 1 is 
an arbitrary small constant independent of L. 

Uniform convergence in the interior [15j: Let {Q ife } be a sequence of global Landau-de Gennes 
minimizers in Aq such that — > as k — > oo. Then (up to a subsequence) Q ifc — > Q° strongly in 
W 1)2 (n, S). 

Let K C be a compact set which does not contain any singularities of Q°. We define 

e L (Q L ) = \\VQ\ 2 + ^. 

Then 

(i) there exists a constant C > independent of L such that 

-Ae L (Q L )(x) <Ce|(Q L )(x) xGK (29) 

for L sufficiently small; 

(ii) we have a uniform bound for e^Q^) in the interior of Q, away from Sq in the limit L — > + i.e. 

e L (Q L )(x) < C'(a 2 ,b 2 ,c 2 ,n) x e K (30) 

for all L sufficiently small and a positive constant C' independent of L; 

(hi) Q ifc converges uniformly to Q° everywhere in the interior of Q, away from Sq. 

lim Q Lfc (x) = Q°(x) uniformly for x £ K. (31) 

A;— >oo 

We emphasize that (|30p and (|31|) only hold in the interior of 0. In Section [5l we extend these 
uniform convergence results up to the boundary for the uniaxial case. 

A uniaxial global Landau-de Gennes minimizer is fully characterized by its scalar order 
parameter s L and distinguished eigenvector n L . The scalar order parameter, s L , is a locally Lipschitz 
function of Q L and hence, is continuous on [27] . From [16], we have that 

< s L (x) < s+ x G tt (32) 

and let 0£ = {x G O; s^(x) = 0} denote the isotropic set of Q L . We have a topologically non-trivial 
boundary condition in ([7]) and hence, every interior extension of Q& must have discontinuities. 
We interpret the defect set of Q L as being the defect set of n L . Let denote the defect set of n L 
and let x n G S^. Then Q L (x n ) = 0, since is well-defined on Q and consequently s L (x n ) = 0. 
We deduce that c8t and from [21], we have that n L is analytic everywhere away from G^. We 
first make an elementary observation about the defect locations as L — > + . 
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Lemma 1 Let So denote the singular set of the limiting harmonic map Q° defined in ([PJ). Let 
x n £ S^. Then 

dist (x n , So) < e(L) 

where e(L) — > as L — > + . 

Proof: Let x n G S 1 ^. As mentioned above, s L (x n ) = and x n G ©£, where @l has been 
defined above. However, the bulk energy density /b(Q l ) converges uniformly to its minimum value, 
everywhere away from So, in the interior and up to the boundary, as L — > 0. Recalling (|28p . we 
deduce that dist(x n , So) — > as L — > + . Lemma[T]now follows. □ 

Lemma [T] is also equivalent to the statement dist(&L, So) — > as L — > + i.e. the isotropic set 
of a uniaxial global Landau-de Gennes energy minimizer converges to the singular set of a limiting 
harmonic map in the limit of vanishing elastic constant. 

Proposition 1 Let Q L be an uniaxial global minimizer of leg i n the admissible space Aq, for 
a fixed L > 0. Then = s L (n L ® n L — |l) for some non-negative function s : — > K + and 
n G Vt /1,2 (0; S 2 ). The following equations hold everywhere in Q, away from the isotropic set 

As L - 3s L |Vn L | 2 = |i (2c 2 (s L ) 2 - b 2 s L - 3a 2 ) (33) 

Anj + |Vn L | 2 n^ + = j, k = 1, 2, 3. (34) 



//ere n^ fc denotes the partial derivative q^-- Alternatively, n L = (sin 6 L cos <p L , s'm6 L sm(j) L , cos 9 L ) , 
where 6 L ,(f> L are functions of spherical polar coordinates (r,9,(j)) centered at the origin. Then L 
and (f> L satisfy the following coupled nonlinear partial differential equations 

V-((s L ) 2 Vd L ) = (s L ) 2 sin# L cosfl L |V/f (35) 
V- ((s L ) 2 sin 2 L V(/> L ) =0. (36) 

Remark: In general, Q G W 1 ' 2 implies that the tensor n <g> n G W 1 ' 2 . However, for simply- 
connected three-dimensional domains, n®n£ W 1,2 (Q) =^ n G W 1,2 (Q; S 2 ) 

Proof: In what follows, we drop the superscript L from Q L for brevity. Since Q is a classical 
solution of (1221). we have that 



Qij,k = d k s ^nmj - + s (nm iife + n^n i)fe ) 

Qy,fcfe = As - ^<%J + 2<9 fc s (nm^fc + n^n^) + s (n,n jitt + n.,n i)fcfe + 2n itk n jtk ) (37) 



where i, j, k = 1 . . . 3, Q^fc = -q^- etc. 
Consider the decoupled equations ([22]) 

1 

3 

and multiply both sides by rij to get the following vector equation 



LQij,kk = t: (2c 2 s 2 - b 2 s - 3a 2 ) Qij 



2 2s 

-n, As + 2d k sn j)k + s (n jykk - |Vn| 2 n.,-) = — (2c 2 s 2 - b 2 s - 3a 2 ) rij. (38) 
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Multiplying both sides of (|38p by n,, we obtain the following scalar equation for s:- 

^As - 2s|Vn| 2 = (2cV - 6 2 s - 3a 2 ) (39) 
3 9-L 



and ([33]) now follows. In ([38]) and (|39|) . we use ([37]) and the relations njrij = l^n^n^ = and 
nin itkk = -|Vn| 2 . 

For ([34"]) . we multiply both sides of the vector equation ([38]) by the derivative nj tP for p = 1, 2, 3 
to get the following system of three equations - 

2d k sn j:P n jtk + sn j>p n j>kk = p = 1, 2, 3. (40) 

Multiplying both sides by the scalar order parameter s, (|40|) simplifies to 

n,- p d fc (s 2 ^- fc ) = p= 1,2,3. (41) 

Next, we note that for a fixed p, (iij, nj tP , ej) form an orthogonal basis at each point x 6 $7, away 
from the isotropic set @l so that 

<9 fc (s 2 n iife ) = \irij + A 2 ej (42) 

where 

Ai = iijdk (s 2 n j:k ) = -s 2 |Vn| 2 . 
We substitute (g2]) into ([38]) to get 



— rij-As - 2s 2 |Vn| 2 rij + A 2 e,,- = (2c 2 s 2 - 6 2 s - 3a 2 ) n, 



from which we deduce that A 2 = 0. Hence 

d k (s 2 n jtk ) + s 2 |Vn| 2 n i = j = 1...3 (43) 

from which (13411 follows. 



An alternative formulation of ([41 p can be obtained by writing the unit-vector field n in terms of 
its spherical angles, 9 L (r, 9, (ft) and 4> L (r, 0, (/>), where (r, 6, (f>) are spherical polar coordinates centered 
at the origin i.e. 

n = (sin 9 L cos cf) L , sin 9 L sin <p L , cos 9 L ) . (44) 



Straightforward computations show that 

= d k 9 L (cos^cos^cos^sin^-sinfl 1 ') +sin6> L d k (\> L (- sinc/> L , cos 4> L , 0) 

d " = d kk 9 L (cos9 L cos<p L ,cos9 L sm(/) L ,-sm9 L ) - (d k 9 L ) 2 (sinfl L cos^ L ,sin6» L sin0 L ,cos(9 L ) + 
+2 cos 9 L d k 9 L d k $ L (- sin0 L , cos 4> L , 0) +sin9 L d kk (j) L (- sin cos L , 0) - sin0 L (9 fc L ) 2 (cos0 L ,sin< 

Substituting ([43]) into ([38]) and taking the dot product of both sides with (cos 9 L cos cos # L sin^ L , — si: 
we obtain 

2d k sd k 9 L + sd kk 9 L - ssin9 L cos9 L \V(j) L \ 2 = 0. (46) 
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We multiply both sides of (|46p by s and equation (J35J) now follows. Similarly, we take the scalar 
product of (|3"5j) with the unit- vector (— sin^, cos cft L , 0) to obtain 

s sin 9 L d kk (ft L + 2s cos 6 L d k 6 L d k (j) L + 2 sin 8 L d k sd k (j) L = 0. (47) 

As above, we multiply both sides of (|47p by ssin# and (|36p then follows. The proof of Proposition [1] 
is now complete. □ 

Comment: We note that for s constant, (34\ ) is equivalent to the harmonic map equations Ano + 
|Vn | 2 n = 0i 

Next, we use asymptotic methods to predict the minimizer profile near isolated isotropic points 
in 6i and establish a 1 — 1 correspondence between isolated isotropic points and isolated point 
defects. 

Proposition 2 Let Q L be an uniaxial global minimizer of leg i n the admissible space Aq, for a 
fixed L > 0. Then Q^(x) = s L (x) (n^(x) <g) n L (x) — Let F n C SI denote the set of isolated 
point defects in n L and let Fl C denote the set of isolated isotropic points of Q L . 

(i) Let £ be an isolated interior isotropic point. Let (r, 9, (ft) denote a local spherical co-ordinate 
system centered at x^;, then 



Vn 1 



a 



as r 







(48) 



where a only depends on 6,(j) and is independent of the radial coordinate r. Then x^ £ r n too. 
(ii) Let x n S T n be an isolated point defect. Then x n G Tl and hence, T n = T^. 



Proof: (i) Consider the coupled equation ([33]): 



As - 3s|Vn| 2 



3L 



2cV - b 2 s - 3a 2 ) 



(49) 



Let X£ S Ti be an isolated isotropic point. Since s L (x) = |Q L (x) (n L (x) <g)n L (x) — ^i) is the 
product of two analytic matrices away from x^, we deduce that s^(x) is analytic for < r < ro, for 
some ro > 0. We are interested in the leading-order behaviour of |Vn L | 2 as r — > 0. 

From the local analyticity of s L , we have the following power series expansion 
s L (x) 



r n g{ 



< r < ri, < ro, n > 1, 



r"g 



o(l) as r 



(50) 
and 



where (r, 9, (ft) is a local spherical coordinate system centered at x^ 

rf, is the radius of convergence of the series (|50p . Further, g ^ for r ^ and {g,h} are analytic 
functions with bounded derivatives in a sufficiently small neighbourhood of x^. Substituting (|50p 
into (|49p and expressing As L in spherical polar coordinates 



As^ 



]_d_ 

r 2 dr 



,ds L 
dr 



+ 



1 







2 C ,L 



we have that 



„n-2 



n(n + l)g + 



r 2 sin 2 i 



+ 



+ 



1 







sin 



cot 



r 2 sin ( 



+ 



ds L 



sin i 



(51) 



2/i r , 

H h /irr + 



r r 2 sin 2 

,i|2 / n 



+ 



+ cot 



3|VnY (r n g(6,<ft) + h(r,9,<f))) 



(r n g + /i) 
3L 



2c 2 (r n o + n) 2 - b 2 (r n g + h) - 3a 2 



as r 



0. 



(52) 
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All the terms on the right-hand side are 0(r n ) whereas the leading order term on the left-hand 
side of (|52|) is 0(r n ~ 2 ). Since h is an analytic function and l^-l = o(l) as r — > 0, we have that 
h r /r, h rr = o{r n ~ 2 ) as r — > 0. Therefore, for (|52p to hold as r — >• 0, we must have 



|Vn L l 2 ~ 



3r 2 



asr^O (53) 



sm' 9 9 

and (USD now follows. It follows that |Vn^| 2 is not defined as r — > and hence, the isolated isotropic 
point Xf, £ T n too. 

(ii) Let x n G T n . Then Q L (x n ) = 0, since Q L is well-defined on ft. Therefore, we must have 
s L (x n ) = and by definition, x n G T^. Combining (i) and (ii), we conclude that = T n . □ 

Comment: By analogy with U8\j . one might expect that uniaxial global minimizers can only account 
for isolated point defects and all higher- dimensional defects are intrinsically biaxial. We hypothesize 
that n = 2 in 150\) i.e. we have a quadratic decay of the scalar order parameter as we approach point 
defects, by analogy with the study of vortices in Ginzburg-Landau theory 

Comment: The estimate fl^#| ) is analogous to a similar result on singularity profiles within the Oseen- 
Frank theory of uniaxial nematic liquid crystals with a constant scalar order parameter s f^. In J2f, 
the authors show that near every singularity x p G ft, we have 

n ~ R^^t (54) 

|X Xp| 

for some rotation R G 50(3). Therefore , 

ir, 12 2 

|Vn| ~ i as x —?• Xp. 

|X Xpl 

The estimate |^<g| ) suggests that we have a similar radial hedgehog-type of profile |5^[ ) near the isolated 
zero xl G ft, for uniaxial global minimizers within the Landau-de Gennes theory. 



5 Far-field results 

In this section, we study the qualitative properties of uniaxial global minimizers {Q 1 '} away from 
the isotropic set @l, in the limit L — > 0. This is equivalent to studying the qualitative properties of 
{Q L } away from the singular set So of the limiting harmonic map Q° defined in ([9]), as L — > + . 

Let Q L = s L (n L ®n i - ^i) be a uniaxial global minimizer for fixed L > 0. Recall that for L 
sufficiently small, 

or equivalently 



< s+ - s L (x) < ei(L) (55) 



Q L l 2 -^4 



< e 2 (L) (56) 
where ei(L),62(L) — > as L — > + , everywhere away from Sq. 



Our first result is an inequality for 



that holds everywhere away from Sq on ft. We do not use Lemma [2] in the subsequent sections but 
keep it as an interesting technical result. 
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Lemma 2 Let A L = ^Q^- fcQfj & &y definition. Then we have the following inequality on Vl\B$(Sq) 
for L sufficiently small 



\A L + \D 2 Q L \ 2 < \ \D 2 Q L \ 2 + a 4 ^, (57) 



where a > 1 is a positive constant independent of L that can be worked out explicitly, Bs(Sq) is a 
small 5 -neighbourhood of So and 5 > is independent of L. 

Proof: The derivation of (J57J) closely follows the methods in @]. In what follows, we drop the 
superscript L for brevity. First, consider the decoupled equations (f22|) ; setting 

f(s) = (2c 2 s 2 - b 2 s - 3a 2 ) 

and differentiating both sides of (j22[) with respect to x p , we obtain 

Qij, k k P = ^ff(s) + /' W^^p ^ P = 1, 2, 3. (58) 

From (|56p and the global upper bound (|23p . we have that |Q| is bounded away from zero on 0,\Bs(So) 
and 

/(*) < o 
f'(s) > 

/"(s) > 0, (59) 
on the set Q \ Bs(Sq), where / (s) = Jj, / (s) = etc. A straightforward computation shows that 

AA = |D 2 Q| 2 + Q ( ,./,Ci, ; . ; , ; , / .. (60) 

and using ([58]) . we obtain 



v^,2 , ,™, 2 /(*) , / W(Q-VQ) 



2 



Ajl = l DQ l + l VQ ll^ + ^ |Q| ' (61) 



From fl59J) and ([22]) . we have the following inequality 



- Ayl + |Z) 2 Q| 2 < |VQ| 2 ^1. (62) 

Finally, we use the inequality 

|AQ| < a \D 2 Q\ 

where a > 1 is a positive constant that can be worked out explicitly. Substituting the above into 

dS2D, 

- A.4 + |D 2 Q| 2 < 2aA 1 -^ < ±\D 2 Q\ 2 + a 4 -^ (63) 
and (15711 now follows. □ 



We recall the uniform convergence result in (|31|) and (|30p . whereby we establish a uniform bound 
for |VQ |, independent of L, everywhere away from So in the interior of Q. The next step is to 
extend this uniform convergence result up to the boundary. To do so, we adapt the small energy 
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regularity theorem in [7] to the Landau-de Gennes framework to obtain a uniform bound for |VQ L | 
independent of L, everywhere away from So up to the boundary. 

Consider a boundary point xo G dCl and define the region O r (xo) = 0(1 B v (xq), where B r (xo) 
is a ball of radius r centered at xo- Let p be a suitably small positive constant such that for any 
xo G dCl, we may choose a coordinate system {x Q } so that xo is at the origin and CI p (xq) corresponds 
to B+(x ) = {x G CI; |x| < p; x 3 > 0}. 

Proposition 3 -Lei {Q Lfe } &e a sequence of uniaxial global minimizers for X eg in the admissible 
space Aq, where Lj, — > as A; — > oo. VFe can extract a subsequence such that Q Lk — > Q° strongly 
in W 1,2 (Cl; S) as k — > oo. Lei xo G 50 6e such that O r (xo) contains no singularity of the limiting 
harmonic map Q°. Then there exists C\ > 0,6*2 > 0, ro > 0,Lq > (aZZ constants independent of 
Lk) so that if 



j 



^|VQ L | 2 + 1 dx<Ci r<mm{r ,p} (64) 

n r (x ) 2 L 

r 2 sup e L (Q L ) < C 2 for all L k < L (65) 
eL(Q L ) = ^|VQ i | 2 + ^p. 



then 



where 



Proof: The first half of the proof of Proposition [3] closely follows the scaling arguments for the 
interior uniform convergence result (|3(jp in [JS] and the second half closely follows the arguments in 
Theorem 2.1 in [7]. 

We first recall from (|26p that since O r (xo) contains no singularity of Q°, 3m(x) G S 2 such that 



Q L (x) - s + ( m(x) ® m(x) - -I 



< e « 1 x G O r (x ) (66) 



for L sufficiently small. 

We continue reasoning similarly to [15]. We fix an arbitrary L/% < Lq. We let < r\ < ^ < 
min ^| and xi G O ri (xo) be such that 

max [% - s) max e Lk (Q Lfc ) = (j ~ n) e Lk {Q Lk ) ( Xl ). (67) 



0<s<^ V 3 / ^( x o) 

Define e\ h = max x£S j ri{xo) e Lk (Q Lk ) = e Lfc (Q Lfe )(xi). Then 



max e Lk (Q Lk ) < 4e^ fe (68) 

^2/3r--ri ( x l) 
2 

where we use the inclusion Cl 2 /3r-r 1 (xi) C 2/3r+ ri (xp), 2 / 3 ^ +ri < ^ by definition of r\ and the 

5 5 

inequalities ([57]) . 



Define r 2 = ri y e^* and let 

R Lfe (x) = Q^ k I xi + I . (69) 
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Let L k = e^ k L k . Then H Lh has the following properties on f2 r2 (0):- 

e Lk (K L ") = ^-e Lk (Q Lk ) (70) 



er. fl 

xefi r2 (o) 



e £fc (R Lfc )<4 e £fc (R L *)(0) = l (71) 



= Y k (2cV " b 2 s - 3a 2 ) R£ (72) 



where s 2 = | |Q Lfc | 2 . 

We next claim that Ti < 1. It is obvious that r2 < 1 implies the conclusion (|65|) . We prove 
this claim by contradiction. Assume that r2 > 1; then using the same arguments as in [7], one is 
led to the existence of a sequence of solutions {R ife } of (|72p on f2i(0) = B^ (0), with the following 
properties: - 

-AB£* + (2c 2 s 2 - 6 2 s - 3a 2 ) r£* = in S+(0) 
3 L k 3 

max e £fe (R Lfc )<4 e £fc (R ifc )(0) = 1 
xesf(o) 



R^ fe |x3=o = Qb I x+ — J== I wii/i 



|VR ifc || X3= o < efc|VQ6| L oo (9 Q), |V 2 R Lfc || X3=0 < ell^Q&U 00 ^) uni/t e k ->• as A; -)• oo 

/ e £ (R ifc ) (fx < o fc -> + as fc -> oo. (73) 
Jb+(o) 

From (j29|) and (fTTj) . we deduce that R Lfc satisfies the following Bochner-type inequality on Bf(0):- 

-Ae Lk (R L «)<C'e Lk (R L k) xeB+(0) (74) 

where C is a constant independent of L k . 

Next, we write R ifc explicitly in terms of its scalar order parameter and leading eigenvector - 



Rf/ = s k (n*nj - Ui^j n k e W 1 ' 2 {Q; S 2 ) 



(75) 



where |R Lfe | 2 = §s| and 



\Sh — S_d < 

' +1 - 10Q 



from (|66p . for A; sufficiently large. From Proposition [H we have that s k and n k satisfy the following 
equations in Bf(0):- 

As k - 3s fc |Vn fc | 2 = -^(2c 2 s 2 - b 2 s k - 3a 2 ) (76) 
6L k 

AnJ + |Vn fc | 2 nJ + 2^1^ = (77) 

| VR i fc |2 = |v Sfc | 2 + 2s||Vn fc | 2 < 4. (78) 
3 
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Prom (|66jl and (|78|). we deduce that 

2Vs k 



|Vn fc | < — . 
s+ 

We combine (|T9"j) and (|77p to deduce that 



Sfe 



< — on B+(0). (79) 



sup |Vn fc | 2 < c5l /3 A; -> oo (80) 
B+ 3 (0) 

where c is a constant independent of fe. In particular, this implies that ||Vn \\L°°(dn) 

< c5l /3 -> as 

A; —7- oo. The proof of (|80p is identical to the proof of Theorem 2.1 in [7J and the details are omitted 
for brevity. 



Next, look at the equation (|33p and introduce the function 

Sfc = S+ — Sfc- 

Then s^ is a solution of the following problem on B^~ (0) :- 

— Asfc = 3sfe Vn 



3L fe 

s fc (x) =0 x G {B+(0) n x 3 = 0} (81) 

where s_ < is a constant. Repeating the same arguments as in [7], we obtain the following 
estimates:- 

Sfe(x) < cix 3 5l /3 x G 5+ /2 (0) 

1/3 

HVsfcH^^n) = ||Vsfc||z,oo( an ) < co^' /e -> oo (82) 
where Co and ci are positive constants independent of L. 

Finally, we define e k = max |o, eL fc (R Lfc ) — (c + co)(5^ 3 |. (Recall that /s(R Lfc ) = on 23 = 
because of the choice of the boundary condition Q&.) Then from (f29j) . ([80]) and ([82]) . we have that 

-Ae fc (x) < C"e fc (x) x G 5+ /2 (0) 

e fe (x) 1^=0 = (83) 
for a constant C" independent of L. Using standard arguments as in [7J, fj83[> implies that 

sup efc(x) < C3<5fc — )• + as k — > 00 

x€B+ 4 (0) 

contradicting (|7ip . The proof of Proposition [3] is now complete. □ 

For the reader's convenience, we quote Lemma 2 from [3] which is used in Proposition 
Lemma 2 from 01/ : Let u{r) be a solution of 

-e 2 Auj + u = on B(0,R) 

to = 1 on dB(0,R). (84) 
Then for e < ¥,u;(r) < e^( r2 ~ R ^ on B(0,R). 
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Proposition 4 Let {Q Lfc } be a sequence of uniaxial global Landau-de Gennes minimizers in the 
admissible space Aq, where — > + as k —> oo. Then as k — >■ oo, we can extract a suitable 
subsequence such that Q Lk — > Q° in C 1,a (0 \ Bs(Sq)) for some < a < 1 and Bs(Sq) is a small 
5 -neighbourhood of the singular set, So, of the limiting harmonic map, Q°, where Q° has been defined 
in (OJ) and 5 is independent of . 



Proof: The proof follows the methods in [3] and the key ingredient is to establish a global bound 
for k , everywhere away from So, for sufficiently small. 

We drop the superscript in what follows for convenience. Consider the equation (|33p on 
Q \ Bs(Sq) and introduce the function 

* = £± ^, (85) 
s+ has been defined in Q and s_ = (^b 2 — \Jb 2 + 24a 2 c 2 ^ /4c 2 . Then (|33|) can be re- written as 

As - 3s|Vn| 2 = -^t/j(s - s_) (86) 
3 

From ([23]) and ((28]), we have that \s\ > |Q| 2 > \s\ - e L where e L as L 0, on Tl \ B S (S ). 
Therefore, 

2c 2 s , , 1 

on fi \ Bs(Sq), where f3 is a positive constant independent of L. 

We note that | VQ| 2 = ||Vs| 2 + 2s 2 |Vn| 2 where Q = s (nfgin — ^i). We recall the global 
uniform bound ([53]) everywhere away from So to deduce that 

|Vn| 2 < C(a 2 ,b 2 ,c 2 ,n) 

onfi \ Bs(So)- Combining the above, we have that ip satisfies the following inequality on O \ Bs(So) 

- (3LAtp + if> < 7|Vn| 2 < D(a 2 , b 2 ,c 2 , 0) (87) 

where 7 and D are positive constants independent of L. Applying standard maximum principle 
arguments, we conclude that 

IMU»(n\B s (So)) <D'(a 2 ,b 2 ,c 2 ,n) (88) 

where D' is a positive constant independent of L. 

Consider the governing equations (|22p for a uniaxial global minimizer Q; they can be written in 
terms of the function ip as shown below - 

AQ = -i- (2cV - b 2 s - 3a 2 ) Q < -atpQ (89) 

oL 



where a > is a constant independent of L, we have used the definition of in (|85p and the uniform 
convergence of bulk energy density everywhere away from So (refer to (|26p). Finally, we combine 
the global upper bound (|23p and the L°°-estimate (|88p to conclude that 

II A QHl°w 4 -(So)) ^ D"(a 2 ,b 2 ,c 2 ,n) (90) 
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where D" is a positive constant independent of L i.e. |AQ| can be bounded independently of Lk 
everywhere away from So. Finally, we use (|90|) and Sobolev estimates to establish {Q ife } — > Q° in 
C 1,a (Q \ B$(So)) as k — > oo, for some < a < 1. The proof of Proposition [4] is now complete. □ 

Comment: One immediate consequence of |23p is that s+ — s < CL, where C is a positive 
constant independent of L, everywhere away from So in the limit L — > 0. This explicitly estimates 
the rate of convergence in 128\) and improves upon a previous estimate in I15f where an analysis 
of the bulk energy density fs in shows that s + — s < C\\f~L where C\ is a positive constant 
independent of L. 

Lemma 3 Let Q L = s (n L ® n — ^i) be a uniaxial global minimizer of leg Aq, for L suffi- 
ciently small. Then for x € O \ B$(So), we have 

\Vs L \ < ei(x) 

||Vn L (x)| 2 -|Vn | 2 |< e2 (x) (91) 
where no and Q° are defined in and ex, e<i — > + as L — > + . 

Proof: Lemma [3] is a direct consequence of Proposition [H Let x G Q \ Bs(So). Then from 
Proposition U we have that 

|Qg(x)-Q9.(x)|<e 3 (x) 

|Qg, fc (x)-Q?- fc (x)| <e 4 (x) (92) 

where Q° is the limiting harmonic map in Q, = and €3,64 << 1. One can directly 

compute 

|VQ°| 2 = 24 |Vn | 2 . (93) 

On the other hand, 

|Q L | 2 = ^ L ) 2 



and therefore, 



3 



q!jQ!jj, = ls L d kS L 



3 

where \s L (x) - s+\ < e 5 (x) << 1 for x G Q \ B 5 (So), from (f28|) . 
One can re- write Q^Qfj- j. as shown below - 

QijQik = (Q§(x) - Q°-(x))Q§ )fe (x) + Q°.(x) (Q§ >fc (x) - Q?- fc (x)) (95) 

since Q?Qy fe = from |Q°| 2 = |s+. Using the inequalities ([92]) . the global bound ([65]) and the 
triangle inequality, we have that 

|Q§(x)Q£. fe (x)|<e 6 (x)<<l (96) 
for x£fl \ Bs(So) and from (fM|) and (|28|) . this necessarily implies that 

|Vs L |<e 7 (L) (97) 

away from So, where £7 — > + as L — > + . 
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On the other hand, from Proposition H] Q L — > Q° in C 1,a (f2; S) as L — > (up to a subsequence) 
everywhere away from So- Therefore, for x £ Q \ Bg(So), 

jlVQ^I 2 - |VQ°| 2 | <e 8 (x) 



(98) 



where es — > + as L — > + . A direct computation shows that 



IVQ 



L\2 



Vs L \ 2 + 2(s^) 2 |Vn 



L\2\n„L\2 



Combining ([28]), ([97]), ® and ([93]), we have that |Vn L | 2 |Vn | 2 as L -)• 0+. Lemma H now 
follows. □ 



Proposition 5 Let Q L = s (n L <g> n L — ^i) 6e a uniaxial global minimizer of X^g in Aq, for L 

sufficiently small. Then for x £ Q\ Bs(Sq), we have that 



s + — s 



9|Vn |' 



L Vo 4 + 24a 2 c 2 



< e B (x) 



(99) 



where eg — > + as L — >• + . 

Proof: Consider the function tp L 



in ([85]) and the equation ([33]) on \ ^(Sq) 



As L - 3s L |Vn 



L|2 



(100) 



where s± = b — ^§ 240-2 c2 ■ Equation (jlOOp can be re-arranged to give 



-LA V 



9|Vn | 2 \ , n 2 s L (s L - S -) 



V6 4 + 24a 2 c 2 



3s L |Vn L | 2 + 



9L 



+ 2c^ 



:A|Vn | 



9|Vn 



V6 4 + 24a 2 c 2 



6c 2 s L (s L 



|Vno| 2 . 



Vfe 4 + 24a 2 c 2 Vb A + 24a 2 c 2 

We note that A |Vno| 2 = O(l) away from So and the right-hand side of (|10ip can be written as 

9L . ,_ ,9 6c 2 s L (s L 



(101) 



3s L |Vn L | 2 + 



V6 4 + 24a 2 c 



:A|Vn r 



= 3^1 Vn'T - 3s+ |Vno| z + 3s+ |Vn |' 
Finally, we use ([28]) and ([9T]) to deduce that 



Vo 4 + 24a 2 c 2 
, 6c 2 s L (s L - 



M |Vn | 2 



(102) 



V6 4 + 24a 2 c 



|Vn | 2 + 0(L). (103) 



3s^|Vn L | 2 - 3s+ |Vn | z < e 
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where eio — > + as L — > + and 



LA ^ 



6c 2 s L (s L - S-) , ,n , ,0 

— ; V y - Vn 2 3s+ Vn 

V6 4 + 24a 2 c 2 ' ' 

-*-) = 
9|Vn | 



as L — > + , since s L — s_ — >• (s + — s_) = Vb 4 + 24a 2 c 2 /2c 2 as L — ?■ + . Combining the above, we 
have that 

' 2 n 9|Vn | 2 



+ /3 V 



Vfr 4 + 24a 2 c 2 



< en 



(104) 



v / 6 4 + 24a 2 c 2 / 

where /? is a positive constant independent of L and en — > + as L — )• + . Proposition now follows 
from the maximum principle and Lemma 2 of [1]. □ 
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6 Generalizations 



This paper focuses on qualitative properties of global minimizers of the Landau-de Gennes energy 
functional on 2D and 3D domains. In the 2D case, we show that the Landau-de Gennes theory is 
equivalent to Ginzburg-Landau theory for superconductors and make predictions about the dimen- 
sion of the defect set, the defect locations and the asymptotic profile of global minimizers close to 
and far away from the defect set. 

In 3D, we focus on uniaxial global minimizers of the Landau-de Gennes energy functional because 
this is the first step in a rigorous study of arbitrary minimizers. The topological defects are contained 
inside the isotropic set of the uniaxial global minimizers. We derive the governing equations for 
the scalar order parameter s L and the leading eigenvector n L ; these equations reflect the coupling 
between the two quantities. We show that the topological defects (or equivalently the isotropic set) 
are necessarily contained in a small neighbourhood of the singular set of a limiting harmonic map 
and establish the vortex-like or radial hedgehog-like profile of isolated point defects. We also study 
the qualitative properties of uniaxial global minimizers away from the isotropic set. In particular, we 
establish the C 1 ' "-convergence of uniaxial global minimizers to a limiting harmonic map, everywhere 
away from the isotropic set, in the limit of vanishing elastic constant. We use this convergence result 
to obtain an expansion for the scalar order parameter s L in terms of L, everywhere away from the 
isotropic set. As mentioned in Section [21 a limiting harmonic map is an energy minimizer within the 
Oseen-Frank theory for uniaxial liquid crystals with constant order parameter. These convergence 
results suggest that Oseen-Frank theory and Landau-de Gennes theory give qualitatively similar 
information away from topological defects and the Landau-de Gennes theory can potentially give 
new information near topological defects. 

As mentioned in Section [IJ some of our results will also apply to uniaxial solutions with bounded 
energy in the limit L — > + . In particular, the strong convergence result in Section U] will hold for 
any sequence of solutions of (|2ip and ()22[) whose energy is bounded from above by the energy of the 
limiting harmonic map Q°, in the limit L — > + . The interior and boundary monotonicty lemmas 
in Section H] hold for all solutions of (|2ip (and hence, of (|22p which is a special case of (|2ip ). In 
particular, (|28p will be valid for all sequences of uniaxial solutions, {Q Lfc } of (|22p . whose energies 
are bounded above by the energy of a limiting harmonic map Q°, in the limit L k ->■ + . We will 
still have the C 1,Q -convergence of {Q Lfe } to Q° as L/% — > 0, everywhere away from the singular set 
of the limiting harmonic map. 

Such sequences of uniaxial solutions with bounded energy do exist, such as radial-hedgehog so- 
lutions on a unit ball with strong radial anchoring conditions [T7]. Radial- hedgehog solutions are 
uniaxial, spherically-symmetric solutions of (|22p . They are analogous to degree +1 vortices in 
Ginzburg-Landau theory and in [17], we use Ginzburg-Landau techniques to study these radial- 
hedgehog solutions, their defect cores and stability properties in the limit L — > + . 

Extensions to biaxial case: Some of the arguments in this paper can be extended to general 
biaxial minimizers of the form ([1]). As an example, let {Q ifc } be a sequence of minimizers (biaxial or 
uniaxial) of the Landau-de Gennes energy functional Teg, in the admissible space _4q. Then {Q Lfe } 
converges strongly to a limiting harmonic map Q° (as in ([9])) in W 1,2 (Q, S) (up to a subsequence) 
[15j . for Lfc — > + as k — > oo. Using the interior and boundary monotonicity lemmas (|24p and (|25p . 
we can show that 

Sb (Q Lfc ) -»• (105) 
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uniformly everywhere away from the singular set, So, of the limiting harmonic map or equivalently 

s-t s+, r -> + (106) 

uniformly away from So, as k — > oo |15j . 

In what follows, we derive the analogue of Lemma [2] in the biaxial case. 

Lemma 4 Let Q L be a global minimizer of T^g , for L sufficiently smaZZp. Let 



Then on Q \ Bs(Sq), we have the following inequality 



1 i o r i2 a A 



L2 



- AA 1 ^ + \ D Z Q L \ <—\D z q, L \ +a 4 — ^ (107) 

Gt mi 

where Bs(Sq) is a small 5 -neighbourhood of Sq and a > 1 is a positive constant independent of L. 



Proof: We start with the relation (|60p 

AA L = |D 2 Q L | 2 + Qfj^Qij^pk 
and drop the superscript L for brevity. 

We need to estimate \Qij t kQij,ppk\ i n terms of |AQ| |VQ| 2 /|Q|. Straightforward but tedious 
calculations show that 

L 2 \Q z3 , k Qi 3 , PP k\ 2 = a 4 |VQ| 4 + c 4 (trQ 2 ) 2 |VQ| 4 + 4b A (Q ip Q P ,, f/ Q^ 9 ) 2 + 
+4a 2 b 2 Q ip Q pj , g Q tj jVQ\ 2 - 46 2 c 2 |Q| 2 |VQ| 2 Q ip Qp i , 9 Q ij - ( , - 2a 2 c 2 |Q| 2 |VQ| 4 + 

+4c 4 (Q • VQ) 4 + 4c 4 (Q • VQ) 2 |Q| 2 |VQ| 2 - 4a 2 c 2 |VQ| 2 (Q • VQ) 2 - 86 2 c 2 (Q • VQ) 2 Qi P Q P j, q Qij, q S 
< C(a 2 ,6 2 ,c 2 )|VQ| 4 (105 

where we have used the Euler-Lagrange equations (|2ip to compute the right-hand side of (|108p and 
the uniform convergence of the bulk energy density to its minimum value away from the singular 
set of the limiting harmonic map. It can be shown that the right-hand side of (|108p vanishes for 
Q £ Qmim where Q m ; n has been defined in ([5]). The details of these calculations are omitted here 
for brevity. 

Secondly, 

L 2^Qj! | A Q|2 = fl 4 |VQ| 4 + 2a 2 6 2 |VQ| 4trQj _ 2a 2 c 2 1 Q | 2 1 VQ| 4 



|Q|2 ' - ' - ' - |Q| 

4 + 3 

27|Q| 



o4 i ^.4 i o 02,^2 Oo3<r> 9ct^ 

-26 2 c 2 trQ 3 |VQ| 4 + c 4 |Q| 4 |VQ| 4 + 26 4 g +r J^, 2 |VQ| 4 > D(a 2 , fr 2 , c 2 )|VQ| 4 



(109) 

where 

D(a 2 ,b 2 ,c 2 ) = 



2 we are not making any assumptions about Q; a general global minimizer for leg in the admissible space A exists 
from the direct methods in the calculus of variations 
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if and only if Q € Q m in- 

Combining ([108]) and ([109]) . we get that 

IVOI 2 

|Qii,fcQy«*| < D'(a 2 , b 2 , c 2 Y-^- I AQ| (110) 

where D' is a positive constant independent of L. Substituting (jllOp into (f60j) and repeating the 
same steps as in Lemma [21 (|107p follows. The proof of Lemma 0] is then complete. □ 

Corollary: Let Q L be a global minimizer of leg (biaxial or uniaxial), in the admissible space Aq, 
for L sufficiently small. Then we have the following interior estimates, away from the singular set, 
So, of the limiting harmonic map Q° in (OJ) :- 

1 | VQ L |2 + < H{a 2^ &2; c 2 ; n) (ln) 

\Q°\-\Q L \<C(a 2 ,b 2 ,c 2 )L on K C £l\B s (S ). (112) 

In particular, the largest positive eigenvalue, \\, of Q L , satisfies the following inequality on the 
interior compact subset K C \ Bs(Sq) 

^-\[<D{a 2 ,b 2 ,c 2 )L (113) 

where s+ /ias 6een defined in (E]) and i/ie positive constants H ,C and D are independent of L. 
Proof: The inequality (jllip is a mere repetition of (|30p : see [15] for a proof. 
Consider the function 

|Q L | = (Qp ? Q^) 1/2 P, ? = 1,2,3. 
Then a direct computation shows that |Q^| satisfies the following partial differential equation 

A | Q l = ^!_(Q_ZQ)! + Q^ r , s = 1 ... 3 (11 4, 

IQI IQI IQI 

where we have dropped the superscript L for brevity. 

On the interior compact subset K C O \ Bs(Sq), we have the following inequalities 

< IQI 2 < \*+ 

|VQ| 2 < Ci(a 2 ,6 2 ,c 2 ) (115) 

where we have used (|23p . (| 105j) and (jllip . Therefore, the first two terms on the right-hand side of 
(|114p can be bounded independently of L. We use the Euler-Lagrange equations (f2Tj) to compute 
the third term on the right-hand side of (|114p i.e. 

Q rs AQ rs 1 , 2|/-v|2 i2, r\3 , 2i/-»|41 IQI f 2\r\\2 ^IQI 2] , ^ 2 |Q| 2 f /^ r Q 3 A 

^q^ = ^l^- g|qi " 6trQ + C >QI ) = -{ c IQI j+vsr^-^ioFj 

We recall from [15] that 

/3 2 (Q) = i-6(^)%[o,i] 
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is the biaxiality parameter and as a direct consequence of (jllip . we have 

/3 2 (Q) = l-6(^) 2 <C 2 (a 2 ,& 2 ,c 2 )L 

on the compact interior subset K C Q\Bg(So), for a positive constant C2 independent of L. Further, 
we have the following sequence of inequalities on K C 0, \ B$(So) 

C 3 (a 2 ,b 2 ,c 2 ) (|Q| - |Q°|) < |c 2 |Q| 2 - - a 2 } < C 4 (a 2 ,b 2 ,c 2 ) (|Q| - |Q°|) 

for positive constants 63,64 independent of L (see (|105|) and (|56p ). 
From the preceding remarks, we deduce that 

A|Q(x)| <a(Q 2 ,6 2 ,c 2 ) + C4(a 2 ,6 2 ,c 2 ) (IQ(x)l ~ IQ0(x)l ^ , xeSl\B s (S ) (116) 
where a is a positive constant independent of L. Define the function 

, (IQ°I-IQI) 

* = L • 

Then using f|116j) . we see that ip satisfies the following inequality on K C Q\Bs(Sq) 

- LA^ + P(a 2 ,b 2 ,c 2 )^j < a'(a 2 ,b 2 ,c 2 ). (117) 

Finally, we apply the maximum principle and Lemma 2 in [3] to deduce that 

\iP(-k)\ <j(a 2 ,b 2 ,c 2 ) xen\B 5 (S ), (118) 

for a positive constant 7 independent of L and (|112p follows. The inequality (|112p improves upon a 
previous estimate in [15] where an analysis of the bulk energy density fg , coupled with (|105|) , shows 
that |Q°| — |Q| < A(a 2 , b 2 , c 2 )\/Z, for a positive constant A(a 2 ,b 2 ,c 2 ) independent of L. 

For (|113|) . we use the following alternative representation formula to ([T]) 

Q L = Sl ®n - + Rl (m 8m-p8p) (119) 

where n, m and p are the orthonormal eigenvectors and 

< s + - S L < C 6 y/Z; R 2 L < C 5 L (120) 

on the interior compact subset K C \ B$(So), for positive constants Cq, C5 independent of L (see 
(|105p and Proposition 7 in [E]). We note that 

IQI 2 = \sl + 2R 2 L 
and hence the inequality (|112p necessarily implies that 

/ od2 \ 1/2 

s + -S L (l + < C 7 L, 
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for a positive constant C-j independent of L. This combined with (|120|) i.e. B? L < C$L yields the 
improved estimate 

< s+ - 5 L (x) < C 8 L x G fi \ B tf (5 ) (121) 

where Cs > is independent of L. Finally, it suffices to note from (|119p that the largest positive 
eigenvalue of Q L is given by 

Af = -S L 

and (fTT3l) directly follows from (fT2Tj) . □ 

One might expect that the techniques in this paper can be generalized to obtain the analogue of 
Proposition [3] and Proposition [5] in the biaxial case. However, this can be accomplished only if we 
have a better understanding of the full Euler-Lagrange equations (|2ip . One strategy is to decompose 
the system (j2~Tj) as follows - 

LAQij = -a 2 Qy - b 2 (Q lk Q kj - ^tr(Q 2 )^) + c 2 Q ij tr(Q 2 ) = 

= (-a 2 - b 2 ^ + c 2 |Q| 2 ) Q y + b 2 O^Qij " QifeQfci + ^IQ| 2 <%) (122) 

where we can think of the first term as being a uniaxial component and the second term as being 
a biaxial component. We need to understand the coupling between the uniaxial and the biaxial 
components and to establish quantitative estimates on the magnitude of the biaxial component, in 
order to derive rigorous results for the structure of global Landau-de Gennes energy minimizers and 
their relation to the limiting harmonic map Q° in ([9]). Other future directions are to characterize 
defects in Landau-de Gennes global minimizers (uniaxial versus biaxial cases), to study qualitative 
properties of Landau-de Gennes minimizers for different choices of the boundary conditions i.e. 
when Q;, ^ Q m in where Q m i n has been defined in © and to study Landau-de Gennes minimizers in 
different temperature regimes. We plan to report on these problems in future work. 
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